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ABSTRACT 


This report presents a technique for calculating the electron transport 
and dose expected behind aluminum shields. The method of analysis was first 
used by Burrell, et al. (NASA SP-169) . This work includes results for a 
source whose angular distribution of flux incident on a shield follows the 
cosine function. Included is a comparison of the dose expected from a cosine 
law source to doses from sources of isotropic and normal angular distribution. 
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ELECTRON TRANSPORT FROM A COSINE LAW SOURCE 


INTRODUCTION 


The purpose of this report is to extend the computational methods of 
Burrell, Wright, and Watts [1) for the calculation of electron transport in 
aluminum shields. The original work considered electron transport for cases 
in which the electron flux was either normally or isotropically incident on a 
shield. It is possible to represent reality more closely be representing the 
flux of incident electrons as following a cosine law distribution. This report 
includes predictions of the el ectron dose rate expected behind space- 
craft shields for each source distribution, normal, isotropic, and cosine law. 


INTERPRETATION OF ELECTRON MONTE CARLO DATA 


To compare this work with previous calculations, it is necessary to 
review the methodology already established. Berger [2] presented electron 
energy transfer data computed by Monte Carlo methods. These data gave 
electron energy and number factors for transmission and reflection at five 
energies between 0. 125 and 2. 0 MeV and at five angles between 0 and 7r/2 as 
a function of reduced shield thickness X = Z/r 0 , where Z is the depth normal 
to the surface and r 0 is the mean range of an electron of energy E 0 . The 
energy transmission and reflection factors are defined by Berger as follows: 


T e (X,0) 


energy transmitted at X 

energy incident on shiel d at angle 0 


U) 


B e (X,0) 


energy reflected at X 

energy incident on shield at angle d 


( 2 ) 


The number transmission and reflection coefficients T.. and B are similarly 

IN N 

defined in terms of the number of electrons transmitted or reflected. 



Particles incident on a shield can penetrate it, be reflected back out of 
it, or be absorbed between the surface and the shield thickness, X. Assuming 
conservation of particles, the sum of the probabilities for the occurrence of 
each >f these events is 1. Thus, the energy or number current is: 

X 

1 = T E tX,0) + B^X.0) + f PtX',0) dX' . (3) 

0 

X 

1 = T n (X,0) + 1 N (X,0) + J N(X',0)dX' . (4) 

0 

The integrals represent the absorption factors to the depth X; PTX',0) and 
N’tX',0) are the probability functions for absorption between depths X' and 
X' •* dX'. Since the effectiveness of shielding is of concern, these probability 
functions for absorption are the quantities for which a description is sought. 

Listing the transmission and reflection factors as a function of X = Z/) 0 
reveals that these factors are very slowly varying functions of incident electron 
energy. For computational purposes, these factors can be considered inde- 
pendent of energy. The energy dependence of transmission and reflection is 
largely taken into account by the energy dependence of r 0 . Because of this 
independence of transmission and reflection coefficients on E 0 , it is reasonable 
to take coefficients averaged over ail source energies in the range of interest. 
The electron absorption factors can then be described as functions of the 
variable X by an empirical curve fitted to Berger's data. 

In a similar manner Burrell, et al. , examined Berger's data for sources 
whose fluxes were normally incident by considering only the factors for 0 = 0 
and isotropically incident by integrating the factors over 0 between 0 and z/2. 
The results for energy and number absorption are listed below, since they will 
be used as a basis of comparison for the cosine law result. 

Normal Incidence 


X 

/ P(X',0)dX' = 1 - [t e (X,0) + B e (X,0)J 


= 0. 95 [l - 


exp (-0. 653X - 2.40X‘ i 


- 6.89X 3 ) J . 


(5) 
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( 6 ) 


/ N(X\0)tlX’ = 1 - [t n (X,0> + B n (X,0)1 
0 L J 

= 0.912 - exp (-0.0512X + 1. 128X 2 - 9. 3bX 3 J . 

Isotropic Incidenc e 

ff/2 X r 

f f P(X\0)dX» sin 0 dO = 1 - It^Xj + B E (X)J 
0 0 

= 0. 72 f" 1 - exp( -2. 47X + 0. 752X 2 

- 11.38X 3 )] . ( 7 ) 

* /2 x r 1 

f f N(X\0)dX» sin d 4t) « 1 - T^X) + B N ^ X > 

A L J 


o o 


= 0. 63 [l - expt -0. 1 52X + 0. 479X 2 

- 11. 84X 3 ) J . (8) 


These expressions yield the desired probability functions for absorption by 
changing variables from X to Z/r 0 and differentiating the equations with 
respect to Z, as follows: 

Normal Incidence 



3 


N(Z,0) 

*0 


/ 

0. 4G7 - 2. 057 | — 

r 0 


+ 


25. 


GO 



exp 


-0. 0512 


/ Z \ 

[ — + 1. 12b 

\ r ° / 



2 


-y. as 



1 10 ) 


Isulro | >H ’ Incidence 




0. 752 
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ANALYSIS OF COSINE LAW SOURCE DATA 

There is physical justification for assuming that a cosine current 
source represents a better approximation to reality than either of the normal or 
isotropic current sources previously examined. The penetration of an elec- 
tron into a shield is a function of its angle of incidence. A spacecraft will be 
randomly oriented to the electron field along the trajectory of the spacecraft . 

It can be assumed that all angles of incidence between the incoming electrons 
and the shielding of the craft occur with equal probability. Thus, the time 


4 


average of the electron flux is isotropic to the surface of the shield. Consider 
a plane-parallel plate shield: the number of electrons that enter a unit area 
cf the shield is proportional to the cosine of the angle between a perpendicular 
to the shield and the vector that represents the incident direction of the 
electron. 

Berger and Seltzer |3] have further reported Monte Carlo energy 

transfer data for such a cosine law source at seven energies between 0.5 and 

f>. t)‘.) MeV as a function of the reduced foil thickness, x. The reflection data 

are limited to the maximum reflection coefficient that is a constant for x 2 0. 3. 

To obtain values of B for x = 0. 1 and 0.2, the original Berger data at 

b 

various angles were weighted by the cosine function, then integrated. The 
cosine law energy reflection coefficients are expressed: 

tt/2 

D 0 (x) ” T f B_,(x') cos 0 sin 0 dfl • (13) 

h, 2 J t 

0 

Berger’s data for number reflection displays considerable variation 
with incident energy E 0 . Despite the fact that some inaccuracy does result by 
averaging over all energies, this method was used to retain the simple 
descriptive nature of the expressions for the number absorption case. This 
procedure reduces tables of data to a more convenient functional form. In the 
number transmission case, however, the evaluation of integrals analogous to 
equation (13) yielded values of B^, so large that, when they were added to 

Berger’s values for T^ and the sum subtracted from one, they gave negative 

values for the absorption coefficients at x = 0. 1 and 0. 2. The creation of 
secondary or knock-on electrons is possible for thin shields of low atomic 
mass. The probability that an electron will cause secondaries increases 
with its incident energy, but the flux of electrons in space diminishes 
exponentially as the electron energies increase. Since aluminum has an 
intermediate atomic mass of 27 and the work is directed toward examination 
of thick shields, the possibility that significant numbers of secondary electrons 
would be encountered was discarded. Using this criteria, the values for 

listed in Table I were chosen. The empirical fit to data using these values 
would appear reasonable. 
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TABLE 1. ENEltGY AND NUMBER REFLECTION COEFFICIENTS 

FI >K C( 161 N E \A W INCTDENC E a 


Vr 0 = x 


U N 

0. 1 

0. 087 

0. 100 

0.2 

0. 110 

0. 150 

0. 3 or larger 

0. 122 

0. 200 


a. Average values foi all energies between 0 and G MeV. 


Equations (14) through (17) were obtained for the energy and number 
absorption factors and probability density functions for a cosine law source. 

Cosine Law Incidence 


”/• X 

/ / P(x», 0) dx’ sin 0 do = 1 

0 0 


T £ (x) + B^U) 


(14) 


= 0.878 1 1 - exp( -2. 08x - 3. 54x 2 - 0. 08x 3 

*/2 X 

f f N(x' f 0)dx' sin 0 d0 = 1 - Tt n (x) + B (x) 1 (15) 

6 6 ^ J 


= 0. 8 [l - exp( -0. 22 5x + 0. 321x 2 - 6. 30x 3 ) . 


pi z, = 4- 
*0 


l. 83 + 0. G22 



+ 1G. 01 




-6. 08 



( 10 ) 


G 




It should be noted that these curve fits are not least squares fits in the 
true sense. Rather, the least squares program devised by Watts and loaned 
to the author fits polynomial coefficients. In this work the data were in the 
form of x, the scaled target thickness versus y, the energy or number 
absorption, i. e. . 


y E (x) = 1 - [T £ (x) + B £ (x)J . 
y N U) = 1 " [ T N U) + B N (> ‘j * 


( 18 ) 

(19) 


lor energy and number absorption, respectively. Thus, in order to obtain a 
fit of the form 


y = A £ 1 - exo ( Bx + Cx 2 + Dx 3 )J 


( 20 ; 


equation (19) is manipulated to give: 


log 


= Bx + Cx2 + Dx3 • 


( 21 ) 


Ii the quantities on the left-hand side of this equation are fitted, the least 
squares program generates the coefficients B, C, and D. The coefficient A 
corresponds to the maximum value for the absorption factor, i. e. , 

A = 1 - (T+ B), (22) 


when T - 0 and B equals its maximum value that is constant for Z/r 0 ^ 0. 3. 
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DOSE RATE 


The absorption probability functions for normal, isotropic, and cosine 
law sources, equations til), til), and 1 1G> , are dependent on Z the shield 
thickness and on E the incident energy through the intermediacy of r 0 the mean 
electron path length. The energy dependence of r 0 is given by Burrell, et al. 

I M »» 


(1.33-0. OlilE) (\l 0. 2713E 2 + 0. 0121 -0.11) 


(33) 


where the factor on the right is an empirical fit to the range data given in the 
tables by Berger and Seltzer 1 4 J . 

The absorption probability function represents the probability that an 
electron of energy F has undergone absorption during passage through a shield 
of thickness Z. An example of the energy dependence of equations (9; , (11) , 
and (1(») is given in Figure 1 for an arbitrarily chosen value of Z = 0. 2 gm/cm 2 . 
The effect of the value of the shield thickness Z on the absorption probability 
for the cosine distribution is shown in Figure 2. 

The energy deposition or dose rate D( Z) obtained from an absorption 
probability function p(Z) for a given differential energy spectrum <p(E) 
incident on a shield is given by: 

10 

D( Z) =K f E<f>(E) p(Z) dE. (24) 

E . 

min 


K is a conversion constant which yields the dose in the units desired. E 

nun 

corresponds to the minimum electron energy necessary to penetrate a shield 
of thickness Z. 

Figure 3 is a comparison of the dose as a function of shield thickness 
predicted by equation (24) for the normal, isotropic, and cosine law absorption 
probability functions. Curves are shown for two examples of idealized 
differential electron spectra. 
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4>(E) 



electrons 
cm*- Me V- see 


( 25 ) 


anil 


<J>( E) = 4. 0 


-4E 

e 


electrons 

cm*-MeV-sec 


(26) 


The graph demonstrates that cosine incident current leads to a larger estimate 
of the dose than that from an isotropic source. This is as expected from the 
comparison of the probability functions in Figure 2. 

The increase in the dose is significant for thick shields to the extent 
that the cosine source is believed to be a better approximation to physical 
conditions in space. If the electron flux is described by equation (25) , the 
dose is as high as a factor of two greater from the cosine source as from an 
isotropic source. 
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MeV ► E (MeV) 


FIGURE 1. ENERGY DEPENDENCE OF THE ABSORPTION PROBABILITY 1 LNCIION FOR 
THREE DIFFERENT SOURCES IN A SHIELD OF ARBITRARY THICKNESS 7 . 



P(Z) 



FIGURE 2. ABSORPTION PROBABILITY FUNCTION OF COSINE SOURCE FOR 

SEVERAL SHIELD THICKNESSES 



DOSE IN RADS 



FIGURE 3. COMPARISON OF ELECTRON DOSE FROM THE THREE 
DIFFERENT SOURCES FOR TWO ELECTRON SPECTRA 
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